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Abstract 



We examine isospin breaking in tiie nucleon wave functions due to the u—d 
quark mass difference and the Coulomb interaction among the quarks, and 
their consequences on the nucleon electroweak form factors in a nonrelativistic 
constituent quark model. The mechanically induced isospin breaking in the 
nucleon wave functions and electroweak form factors are exactly evaluated in 
this model. We calculate the electromagnetically induced isospin admixtures 
by using first-order perturbation theory, including the lowest-lying resonance 
with nucleon quantum numbers but isospin 3/2. We find a small (< 1%), but 
finite correction to the anomalous magnetic moments of the nucleon stemming 
almost entirely from the quark mass difference, while the static nucleon axial 
coupling remains uncorrected. Corrections of the same order of magnitude 
appear in charge, magnetic, and axial radii of the nucleon. The correction 
to the charge radius in this model is primarily isoscalar, and may be of some 
significance for the extraction of the strangeness radius from e.g. elastic for- 
ward angle parity violating electron-proton asymmetries, or elastic '^He{e, e') 
experiments. 

PACS numbers: 13.40.-f, 13.40.Ks, 14.20.Dh 
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I. INTRODUCTION 



Parity-violating electroweak lepton-nucleon scattering such as N[e, e')N, N{v, v')N is 
of particular interest to subnucleon physics because it is sensitive to the Standard Model 
fundamental coupling constants a, sin^ 6w and various polar- and axial- vector current form 
factors of the nucleon [|l|. The standard analysis, that led to the current interest in these 
processes, is based on the following assumptions: (i) Lorentz + translational invariance, 
(ii) Standard Model of lepton and quark electroweak interactions, (iii) one boson exchange 
approximation, (iv) significance of only u and d quarks in the nucleon states, (v) good 
parity of the nucleon states, and (vi) good isospin of the nucleon states. In view of pending 
A^(e, e')N experiments 0, and of pending and completed elastic N{v, h'')N measurements, 
it seems worthwhile to explore the corrections due to the relaxation of some, or all, of these 
premises. So far, assumptions (ii-v) have been relaxed and their consequences explored 
|T|,^. In this note we will relax assumption (vi) in a specific quark model of the nucleon and 
look at its effects on electroweak form factors of the nucleon. The results are of immediate 
relevance to the interpretation of the isoscalar axial coupling as measured e.g. in elastic 
neutrino- nucleon scattering They will also be relevant for isoscalar vector couplings 

extracted in future neutrino and parity violation experiments. This is because intrinsic 
isospin breaking, although expected to be small, modifies observables in much the same way 
as nonzero strangeness content of the nucleon does. 

A common present-day picture of the nucleon is one of the ground state of three 
quarks bound by gluon exchange according to quantum chromodynamics (QCD). Given the 
intractable nature of this strong coupling few-body problem, we resort to a nonrelativis- 
tic constituent quark model with harmonic oscillator confining quark-quark interactions, 
chosen for its simplicity. The isospin breaking (IB) in our model is due to the following 
two sources: quark mass differences, and electroweak interactions, which is essentially the 
same as in QCD Since both of these corrections are expected to be small relative 

to the isospin-conserving components, we will evaluate each one in its own right without 
considering the cross-terms. The results of our analysis will justify this assumption ex post 
facto. Furthermore, we will confine ourselves to the investigation of electroweak interactions 
among quarks in the nucleon and neglect their self-interactions. The results will turn out 
to be particularly sensitive to the shift in the nucleon's center of mass, which is faithfully 
reflected in this model. 

In this paper, we begin (Section II) with a general discussion of isospin breaking and the 
connection to electroweak form factors. The purpose is to formalize and clarify the extent to 
which e.g. the strangeness content of the nucleon can be distinguished from IB modifications 
in a model independent way. In section III, we describe the specific quark model we use, and 
outline the calculation of IB effects due to electromagnetic quark-quark interactions, and 
quark mass differences. In section IV, we present the results of these calculations and discuss 
them. We find that some IB effects are significantly larger than one might naively expect. 
Only one of the corrections (the weak neutral current charge radius) we have found appears 
to present possible qualitative problems for extracting non-trivial strangeness content, if that 
should turn out to be present. The calculations presented in sections III and IV are model 
dependent, and are manifestly just a part of a larger set of corrections. We conclude section 
IV with a brief discussion of future directions to be pursued. In section V, we summarize 
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our results and draw our conclusions. 



II. MODEL-INDEPENDENT FORMALISM 

One can define isospin at two different levels: (a) at the level of quarks, where the 
Standard Model is constructed; and (b) at the hadron (here, nucleon) level, where all of 
experimental data originate. We will need both levels to define IB in our formalism. That 
fact limits the applicability of the present formalism to gnar/c models of the nucleon and leaves 
out nucleon models with purely hadronic degrees of freedom such as Skyrme's. Attempts to 
define isospin-breaking admixtures to the nucleon wave functions exclusively at level (b) can 
easily lead to tautologies and circular arguments. In order to avoid that we will first review 
the consequences of exact isospin, and then introduce notation which should clearly denote 
the origin and character of the various contributions to the total (physical) form factors. 

In the isospin-symmetric limit the isospin transformation properties of the operators are 
identical to those of the matrix elements, as a direct consequence of the Wigner-Eckart 
theorem. The main consequence of isospin breaking in the nucleon wave functions is the 
loss of these, previously exact, isospin transformation properties of the electroweak current 
matrix elements as compared with those of the current operators themselves. A simple 
example will best clarify this statement. The third component of the isovector axial current 
in the "nuclear domain", where only u, d quarks are kept, reads [| 

JU = ^Y^/^Tsg • (1) 

When sandwiched between two nucleon states |A^) with "impure" isospin, this does not yield 
just +\Fa for the proton and —\Fa for the neutron, but rather it induces a small but finite 
effective isoscalar axial form factor. To see the relevance of these comments, we turn to the 
weak neutral current (NC) J7^*" = — ■ The axial neutral current, in the presence of 
three flavors, is 

= \ (w7/.75M - d'lf.'lbU - S7^75s) . (2) 

consists of the nuclear domain isovector current and a strange quark part. It is these 
"strange" (hidden strangeness- induced) terms that are a topic of considerable recent interest. 
They would also yield an effective isoscalar axial form factor. The hidden strangeness and 
intrinsic isospin breaking effects cannot be separated by experiment, thus any nontrivial 
isospin breaking could potentially affect a determination of the strange form factors. 

Similar, but not identical comments hold for the charge and magnetic form factors, 
because both the isoscalar and isovector current operators enter the EM current. In the 
following we will define the effects of such isospin-breaking in the nucleon wave functions on 



^Note that our conventions on the definitions of electroweak currents differ slightly from ref. 
We adopt the Bjorken and Drell (BD) metric and conventions, with the exception of the 
normalization of the Pauli form factor which we take to be i*2(0) = k. 
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electroweak form factors. We define the axial- vector current form factors for real nucleon 
states as follows: 



(iV(p')b(M7/.75« - h^-i^d)\N{p)) = u{p') 



/u—d 



u{p) 



(3) 



where ± = {N\t'%\N) for proton or neutron matrix elements, respectively, and — Qq — 
is the invariant momentum transfer. The notation for the form factor superscripts is as 
follows: the upper left-hand index refers to the quark operator on the left-hand side side 
of the equation. The upper right-hand index refers to the nucleon isospin operator. In 
the above case, the quark operator is u — d and hence quark-isovector. It induces a large 
nucleon-isovector '^-^pp-'" axial form factor and a small nucleon-isoscalar axial form factor 
u-djpp+n correction. We also define the "strange" form factors Si^2,A as follows 

{N{p')\is-f^75s)\N{p)) = u{p') \SA{q'h,7,] u{p), 



{N(p')\(s^,s)\N(p)) = u(p') S,{q')^^ + 



iCTnuq 
2M, 



N 



S2{q') 



u{p) 



(4) 
(5) 



which are all isoscalar to first approximation. We ignore here any (isovector) IB corrections 
to the strange form factors of the nucleon, since the strange form factors themselves can 
be viewed as a consequence of flavor SU(3) symmetry breaking admixtures and hence are 
expected to be small relative to the u, d induced ones, thus making the IB corrections to 
them "doubly small" i.e. a second order effect. 5'i(0) = follows from the fact that the 
nucleon has zero total strangeness. 

With our conventions, the EM current operator reads 



tEM 



^{u^u.u - d-fi^d) + ^{u-fi^u + d-f^d - 2s-f^s) 



The polar- vector part of the weak neutral current operator is 



J. 



■NC 



1 



{u^^u - d^^d) - -s^i,s - 2 sin 



EM 



(6) 



(7) 



where 9w is the weak mixing angle. We define polar-vector form factors for real nucleon 
states as follows: 



{N{p')\-{u^^u-d^,d)\N{p)) = 



N 



{N{p')\-{uj^u+dj,d)\N{p)) = 



u{p') 



ydpp+n u+d pp-n-^^ _|_ }_^+dpP+n _|_ u+dpp-n-^ fipq 

^ 2 2Afjv 



u{p) 



u{p) 



(8) 



(9) 



In this notation, if isospin were a good symmetry of the nucleon states, and there were no 

strangeness content, then _ u+dpp-n^q2^ _ g ^j^g Wigncr-Eckart theorem), 

and "+<^i?p+"(^g2j -y^^Q^i^j ^j^g usual nucleon isoscalar form factor, -y^^Q^jj 
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the usual nucleon isovector form factor.^ With the above definitions, we can write weak 
matrix elements in terms of the above electromagnetic and the strange ones with some, as 
yet unknown, nucleon IB corrections. For example, 

F^^ = ^r^Fr"±"-^Fr")-2sin2^H^i^^/-l^, , (10) 

where i = 1,2, i.e. for polar- vector neutral current form factors. Similarly, we have for the 
axial- vector NC form factors 

F^% = ^("-'^Fr- ± "-'^Fr") - ■ (11) 

Taking the difference between the proton and the neutron polar-vector NC form factors, we 
find 

- i^T = (1 - 2 sin^ Ow) (i^™ - i^T) - "^'^r" • (12) 

Here the isospin breaking form factor '^+'^p^~^ appears as a correction to the usual isovector 
form factor relation, and strangeness content cancels out. In principle, this means that the 
IB isovector correction '^+<^p'P~'^ is measurable to this order in perturbation theory, given a 
complete set of experiments, i.e. if we had both the proton and the neutron experiments. 
Similarly, we find for the isoscalar part 

Fl"" + Ff^^ = -2 sin^ e^{F^^ + Ff^) - S, + ^-^F^ ■ (13) 

The final term again represents isospin breaking effects. For the axial form factors we find 

FS - ^An = "-'i^r" , (14) 

FlS + FZ = "--'FT - Sa . (15) 

As discussed above, the isospin correction terms in Eqs. (p!3| ) and ( [I5| ) enter exactly like the 
strange form factors do. All of these IB corrections are expected to be extremely small, of 
course, and will be estimated in this paper. 

Finally, there are some constraints on the IB admixture-induced form factors at zero 
momentum transfer stemming from exact symmetries. Gauge invariance says that the elec- 
tric charge of the proton must be unity, so -Fi^p(O) = 1, and Fi^„(0) = 0. In fact, both 
«+d^p-n.j^g-j "~'^Ff^"(0) vanish identically at zero momentum transfer to leading order 



in IB interactions; the first relation is due to the Ademollo-Gatto [IT2| theorem, the second 



due to baryon number conservation. The latter result follows from the orthogonality of 
excited states used in the first-order perturbation theory and the fact that the charge form 
factor reduces to the norm of the ground state, at zero momentum transfer. This statement 
does not apply to the magnetic or axial form factors at any g^, nor to the electric one at 
nonvanishing momentum transfers. Such correction terms will be typically of order a. In 



^Note that we have absorbed a factor of 1/2 into the definition of the quark isovector, "u-d", 
superscript, and a factor of 1/6 into the quark isoscalar, "u+d", notation. 
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the next section, we calculate, in a simple constituent quark model, the IB corrections to the 
Walecka-Sachs ^ form factors ''+'^0^'' (q^) and i = E,M, which are related 

to the Dirac and Pauli form factors F12 by 



GM{q') 



F,{q') + F,{q' 



(16) 



III. THE MODEL 

We use the constituent quark model []^,|^ to calculate both of the above mentioned 
isospin-breaking corrections and use the harmonic oscillator Hamiltonian for that purpose. 
This three-body problem can be reduced to two uncoupled harmonic oscillators by applica- 
tion of the (equal mass) Jacobi three-body coordinates 

P = ^ (ri - rs) (17a) 

A = -^(ri + r2-2r3) (17b) 

R = ^ (ri + r2 + ra) . (17c) 

Then, the Hamiltonian consists of two independent harmonic oscillators with equal mass m 
and a freely moving center of momentum with mass M = 3m, 

Pi 3fc / 2 2\ _ (18) 
2M 2m 2m 2 ^ ^ ^ 

Solutions to the Schrodinger equation with this Hamiltonian are well known and have been 
tabulated in |^,§,|1^ for low-lying nucleon resonances. A simple confining interaction such 



as the harmonic oscillator leads to spatial wave functions that are Gaussians, in the case 
of the ground state, or that decay as Gaussians at large distances, for any other state of 
the system. That, in turn, leads to Gaussian EM and axial form factors which fall off far 
too rapidly at high values of momentum transfer compared with experiment. This model 
is clearly rather simple, but should be adequate for the purposes of identifying qualitative 
features, and making first estimates of small IB effects. In this spirit we have neglected the 



strong-hyperfine-interaction interference with the IB terms in the Hamiltonian [[10| , |14| . The 



former is an important part of the extended constituent quark model [0, and its effects 
have been evaluated for simple IB observables such as the hadron mass differences [lT0| , p!4 . 
But, in our case it complicates the evaluation of the observables to such an extent that we 
relegate its inclusion to the future. 



A. Mechanical Corrections 

The notion of quark mass is an ill-defined one: free quarks have not been observed. 
Traditionally one distinguishes between two kinds of quark masses: (1) the current quark 
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mass TUq which is one that the quarks would have in the absence of all strong interactions; 
(2) the constituent quark mass mg which is the measure of inertia of a quark moving within 
hadrons. There is of course considerable difficulty in precisely pinning down the values 



of the current quark masses [0, and even more in connecting these to constituent quark 
masses. Weinberg has used chiral symmetry and empirical meson masses to argue that 
rriu — 5MeV;mrf ~ 9MeV;ms ~ 140MeV and mq = rriq + const, where const ^ 330MeV. 
This, in turn, implies 

Anig = rriu — rrid = Amq = mjj — m£, ~ — 4MeV, (19) 

and those are the values that we useQ Consequently, we expect the leading isospin breaking 
corrections due to explicit quark mass differences (which we call "mechanical corrections") 
to be of 0{AmQ/mQ), where ArriQ/mQ ~ —1/85 i.e. rather small. They can be exactly 
evaluated if one assumes harmonic oscillator confining quark-quark interactions with unequal 
masses, 

3 2 u 3 

where m = rrti = m2 7^ = m . The first two (p and A) of the three-body Jacobi 
coordinates in this case are the same as Eqs. ([T7a]) and ([T7b| ), but now carry quark number 
indices e.g. 

P3 = -7^ (ri - 1-2) (21a) 



V2 
1 

7q 



(ri + r2-2r3) , (21b) 

o 

and the center of mass is shifted to 

R = — — ^ (mri + mr2 + m'rg) . (22) 

[Zm + m j ^ ' 

This allows a separation of variables and an exact solution to the problem in terms of two 
harmonic oscillator wave functions with two different masses, 

"iram! , . 

2m + m' 

and a freely moving center of momentum with mass Mjv = 2m+m' (where m = ruu, m' = nid 
in the proton, m = rrid, m' = rriu in the neutron): 

H = ^ + ^ + J^ + ^(p' + X^). (24) 
2Mn 2mp 2mx 2 y ^ ^ ' 



3t 



Isgur [10 1 finds Am^ = — 6MeV in his extended version of the present model that includes the 
strong hyperfine interaction effects. 



7 



The total wave function factors into the CM plane wave solution and an internal motion 
wave function 



^^(P,)) = (^)'^'exp («(P, ■ R - E,t)) |iV(940)) 



where E}^ = + | {up + ux), lu = -y/S/c/m, and the factor 



(25) 

is the square root of the 

inverse Jacobian for the above three-body Jacobi coordinates. The spring constant k is an 
adjustable parameter of this model. 

For spin- and isospin-independent quark potentials the complete internal wave function 
factorizes into a product of the spin, isospin, and spatial wave functions. Once the isospin- 
dependent terms are introduced into the Hamiltonian, however, the factorization property 
breaks down and the spatial and isospin parts of the internal wave function become entangled 



|A^(940)) 



The spin parts are given in a standard notation: 



(26) 



XT 



V2 
1 



(a/3 — Pa) a 



X| = {2aa/3 — a(3a — (3aa) , 
while the spatial-isospin parts for the proton are given by 



= 

p 


{i)a{udu) - 


- tpQ^duu)) 




^p = 

n 


{i)Q{dud) - 


- ipoiudd)) 




$^ = 

P 


{2ipo{uud) 


— ipQiudu) - 


- ipo{duu)) 


<|)^ = 

n 


{2il)o{ddu) 


— tpoldud) - 


- ipoiudd)) 



where 



TpQiudu) 



opUJp 
7T 



\ 71 / 



exp 



1 /p| A|- 

2 Wp Rl. 



mu^JuV^"^ frUdUJdV/'^ /m„,a;,■^^/^ 



TT 



\ TX J 



udu 



® udu 



(27) 



(28a) 



(29) 



X exp ( -- (^u;„m„(ri - R) + ujd'md{r2 - R) + u;„m„(r3 



R) 



the second line of Eq. (g^) is meant in Faiman and Hendry's |]18[ symbolic sense only and 
is not to be used in actual calculations. Here = rUaUJa = y/Sknia, a = p, X. The 
corresponding neutron expression is found by switching u with d everywhere. The simpler 
isospin-s?/mmetnc wave function can be read off directly by setting = = m. 



8 



It factorizes into three parts as expected. The complete wavefunction is crucial to our 
subsequent work, however, as straightforward use of the standard factorized nucleon wave 
functions fails to correctly treat isospin breaking. 

We are primarily interested in the elastic charge and magnetic form factors Ge{(1^), 
GmIq^) defined as the Fourier transforms of the charge and current densities 



dR(^iv(P')|JEM(R)| ^Jv(P)) exp(zq • R) 



(27r)35(P' -q-P)e 



(30) 



in the nonrelativistic limit, where I'^n) is the ground state of the exact Hamiltonian Eq. (p4D, 
P and P are the CM momenta of the initial and final state nucleons, respectively, and 
q = P - P is the three-momentum transfer. We work in the Breit frame defined by 
q = 2P = — 2P, which ensures that E = E . In the Eq. (BDI) above, R is the photon 



position vector and not the CM position vector, S and 1 are the Pauli and unit matrices, 
respectively, operating in the nucleon spin space and (S)Ar, {1)n are their matrix elements 
taken between nucleon spinors. The quark EM current density operator is the nonrelativistic 
reduction of the Dirac fermion current 

Jem(R) = E ^ (-^ {Vr., 5(r, - R)} + X [Vr, 6{r, - R)]) , (31) 

i=l * 

where = |(| + r^) is the quark electric charge operator. Similarly, for the axial current 
in the isospin-symmetric limit, we have 

rfR(^^(P')|Jl(R)| ^^(P))exp(zq-R) 

= {2nf6{P' - q - P)(S)^ (ir^)^FA(q2) . (32) 

Here, tn are the isospin matrices for the nucleons, and the quark axial current reads 

3l{K) = j:^^^S{r,-K) . (33) 

i=i ^ 

Thus we have made yet another assumption: our constituent quark electroweak interactions 
are identical to those of the corresponding current quarks. Specifically, this means that we 
assume vanishing anomalous magnetic moments and no form factors for constituent quarks. 
For a discussion of this point see section IV. C. 

Despite the exact solvability of the Schrodinger equation with the Hamiltonian Eq. (^4|), 
we will expand the exact form factors in powers of q^. The reason for this is that a simple 
confining interaction such as the harmonic oscillator leads to Gaussian spatial wave functions, 
which in turn lead to an unrealistically rapid form factor fall-off at high q^. Therefore, it 
only makes sense to talk about the leading terms in the expansion in powers of momentum 
transfer. 



9 



B. Electromagnetic Corrections 



The electroweak interactions among quarks are readily divided into parity-conserving 
electromagnetic (EM) terms that are of 0{a), where a ~ 1/137 is the fine-structure 
constant, and the parity-violating weak interactions that are of 0{Gf) where Gp — 
(1.023 ± 0.002) X 10~^Mp is the Fermi weak coupling constant. The latter terms lead 
to parity-admixture corrections that have already been analyzed in Ref. IQ. Since we are 
working in the nonrelativistic approximation, we expand the invariant M0ller operator in 
powers of momentum over mass, which leads to the Coulomb interaction as the leading EM 
term 

^EM=E ■ | . (34) 

The hyperfine and the spin-orbit coupling interactions are "higher-order" corrections in 
— . The Coulomb interaction Vem induces a small 0{a) abnormal isospin admixture in 
the nuclear wave function. To determine this admixture we use the first-order (Rayleigh- 
Schrodinger) perturbation theory in the perturbing Coulomb potential, as applied to the 
Hamiltonian 

H = Ho + Vem . (35) 
The ground state of the nucleon l^^o), to 0{a), is given by 

|$„) = |*„) + j:i*,.)<MlillM + 0(«^), (36) 

where are the exact eigenstates of the isospin-symmetric Hamiltonian Hq: Ho\^n) = 
En\^n)- This involves evaluating the isospin-breaking admixture to the nucleon wave func- 
tion defined by the admixture coefficients e„: 

= (MMM (37) 

The sum extends over alln, i.e., over infinitely many excited states of the nucleon. 

The isospin-admixtures in the wave function generate isospin breaking contributions to 
the elastic nucleon electroweak current matrix elements. These are determined by using 
the first-order perturbation theory parameters e„ and the vector and axial current tran- 
sition matrix elements between the nucleon and its excited states calculated in the non- 
relativistic impulse approximation. The definition of the isospin-violating corrections to the 
non-relativistic current matrix elements to 0{a) is thus 

(^o|J| ^o) = ('^'olJl'^'o) + 

E -ET^-ET ((<fo|J|$n)($n|V^EM|$o) + ($o|'^l5M|$n) ($n| J |$o)) , (38) 



Eq — En 



and similarly for the charge density elastic matrix element (^&o|p| ^o)- These formulae have a 
simple Feynman diagrammatic interpretation shown in Fig. |l|. It is important to remember 
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that the above 0{a) corrections do not constitute a complete set. There are other graphs, 
e.g. those shown in Fig. ^ that contribute to the same order in a, but are not included 
in Eq. |38|. Our class of corrections, however, is gauge-invariant and therefore a physically 
sensible subset. 



IV. RESULTS 

The first and foremost question is: do we observe any change in the static properties 
of the nucleon, i.e., in the zero momentum transfer values of the form factors? Then, the 
second question is: how do the effective radii of the nucleon change? Both questions will be 
answered separately in each case. 



A. Mechanical Corrections 

1. Corrections to static nucleon moments 

It turns out that two, the charge and the axial, out of three static couplings remain 
unchanged. In the case of the isoscalar "electric" form factor (convection part of the vector 
current) at zero momentum transfer this is a straightforward consequence of baryon number 
conservation. The isovector electric form factor at zero momentum transfer is unrenor- 
malized as well, in agreement with the AdemoUo-Gatto theorem The axial coupling 

constant also remains unchanged, but there does not seem to be as deep a reason for that as 
for the charge conservation. In other words, this quark model as it stands, does not predict 
the existence of an additional isovector or any isoscalar axial couplings due to the Coulomb 
interaction between the quarks and/or different quark masses. The magnetic moments, 
however, are corrected. 

The change of variables from the equal mass three-body Jacobi coordinates Eqs. ( |17a| - 
c) to the unequal mass case (Eq. amounts to a shift of the center of mass (CM) of 
the system. That shift provides the main source of mechanical corrections. The change of 
the oscillator frequencies plays a secondary role (see tables I and II). E.g., the anomalous 
magnetic moments are renormalized due to the change in the quark and nucleon masses, via 
the definition of the magnetic moments: 

{N T I^^J^SliV T) = E(iV T \^CT,\N T), (39) 

where eGM,Ar(0) = + ck^. It is important to note that the quark mass mj is now an 
operator in isospin space mj = m + Am|r^. Then we find 

^ ^ ^ ^ 8mrf + niu ^ 1 A ^ '^rnu\ A ^ 8md\ ^ 
^ 3mA,p 9 V nid J \ J 

2\ md + 2mu 2 / 2md\ ( 2mu\ , , s 

3/ mA,n 9 V vnu J \ vrid J 

Corrections to nucleon static electroweak moments are summarized in Table I. Since the 
magnetic moments are the only static nucleon properties that do receive corrections in this 
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model, we will separate them into various parts according to the nomenclature defined in 
Sec. II: 



Gm(0)p+" = 


1 

9 


7 + 


rriu 
— + 


rrid' 
rriu. 


~ 1.0 


Gm(0)^-" = 




3 + 


ruu ^ 
rrid 


rrid 
rriu 


~ 5.0 


Gm(0)p+" = 


1 
3 


'rrid 
mu 


_ rUu 
md. 


~ 0.008 


G'm(0)p-'^ = 


1 
3 


'rrid 
mu 


_ rUu 
md. 


~ 0.008. 



(41a) 
(41b) 
(41c) 
(41d) 



We can immediately use the results from Eqs. ([T^ T3) to find 
1 



M,n 



(0) 



2 
1 

2 

1 r 

2 - 
1 r 



1 - 4sin2^H/)G™ (0) - (0) - 5m(0) - "+'^Gm(0)^-" + "^'^Gm(0)p+" 

(42a) 



(1 - 4sin2^H.)G™ (0) - (0) - 5m(0) 



1 - 4 sin^ ew)GlfM - Gl^M - 5m(0) + "-^Gm(0)p+'^ + "+'^Gm(0)^-" 



(1 - Asin' dw)Gl^M - Gf^;,(0) - 5m(0) + 2"-'^Gm(0)^+" 



(42b) 



where the second lines of these two formulas represent this model's results. Using the values 
in Eq. ( [41c[ -d), the proton weak magnetic form factor is essentially unaffected, while the 
neutron weak magnetic form factor gets a small but nonzero isospin-breaking correction of 
approximately 0.3%. We conclude that, in this model, the proton magnetic moment mea- 
surements in elastic parity-violating electron-nucleon scattering can he directly interpreted 
as nucleon strangeness content, whereas the neutron ones are subject to the above small IB 
correction. It is not clear, however, to what degree the above cancellation of IB terms in 
G^Pp is model dependent. 



2. Corrections to nucleon radii 

The dependent isospin-breaking effects are present in all of the form factors. A few 
words about their evaluation are in order. The nucleon wave function does not factor into 
three coherent parts due to isospin breaking: the isospin and the spatial wave functions are 
now coupled (see Eq. |28a| ) and that has to be taken into account. With proper book-keeping 
of quark flavors in the spatial wave functions, one derives the following formulae: 



2M„ 



eG"" 

qEM 



'M,p 



2Mr 



2ed{d)n + eu{u)n 



— {u)p + —{d)p 
rUu rud 

1 2 

(m)„, H {d)n 

niu nid 



(43a) 
(43b) 

(43c) 
(43d) 
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-I md)n + {U)n] 



(43e) 
(43f) 



Here {u)p is defined as the Fourier transform of the spatial wave function matrix element of 
a u quark in the proton i.e. {u)p = (expiq ■ r„)p, etc.. They are evaluated as 



{u)n 
{d)n 



(expiq-r„)p 

(exp iq ■ rd)p 
(exp iq ■ r„)„ 
(exp iq ■ Yd)n 



exp 



exp 



exp 



exp 



Rpp + 3 




_3q^ 
2 

_3q^ 
2 



Mp 

Mr,. 



R 



Xp 



R-Xn 



(44a) 

(44b) 

(44c) 
(44d) 



Upon inserting these into Eqs. ( |43a| -f) and expanding in Taylor series, we find the results 
that are shown in Table 2. They all turn out to be small (< 1%), and hence are not likely 
to be distinguishable from the experimental uncertainties in the "known" EM form factors 
form factors. 

One case is of particular importance, though, because of the vanishing of the strange 
and neutron electric form factors in the static limit: the weak NC charge radius. Since 
the nucleon strange radius is then the first term in the Taylor expansion of the nucleon 
strange form factor, it is of immediate experimental interest and will be measured in e.g. 
parity-violating (PV) elastic ^He{e, e') and forward angle p{e, e')p experiments. It is hence 
important that we know the relevant IB corrections. As with the magnetic moment in the 
previous subsection, we break up the total correction to the nucleon EM charge radius into 
its isoscalar and isovector components 



u—d /2 \ n 



m 



1 



2 \ M 



9 / mu \ 
"2 \Mp) 

3 /mrfV 
2 \mJ 



' 3 



RL + 3 ^ 



P2 

^Xn 



+ T RL + 3 



o ITT ^Xn - T 



2 \MnJ 



4 



(In the isospin symmetric limit, e.g., (r'^y 
numerical values into these formulae gives 




Rpn + 3 



(45a) 

(45b) 

(45c) 
(45d) 

R^.) Inserting the 
(46a) 
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,_rf^^|^p+„^^2 ^ 2.1% (46c) 
^ 0.1% . (46d) 

There is a simple physical picture which roughly explains these numbers. The bulk of 
the effect arises from the shift in the position of the center of mass when the constituent 
quarks are given different masses. In the case of the EM charge radius, we are evaluating a 
weighted sum of the quark positions squared. The isovector rms charge radius is weighted 
by the sign of the quark charges, and thus in the neutron, although the shift of the center 
of mass is in the opposite direction from that of the proton, the sign of the quark charges is 
also switched, resulting in the same increase in radius for proton and neutron. For the same 
reason the isoscalar rms charge radius is very small. Thus the bulk effect is primarily due 
to ^-<i _ This simple geometrical picture yields A(r)^/i?^ ^ 2{AmQ/3mQ) ^ .8% for 
proton and neutron separately, and these are then additive. The remainder arises primarily 
from the modification of oscillator frequencies due to reduced mass effects, see Eq. (|23|) . We 
see that the complete IB in the rms charge radius is due to fairly simple kinematical effects, 
which, however, can only be evaluated if the CM motion is correctly accounted for. 

When we include strange form factors, we find 

(ri)?^ = I [(1 - 4sin2 ^iy)(ri)™ - (ri)™ - (r|). - -^'^{rly-- + -'^(r|)^+"] , (47a) 

(4)^ = I [(1 - 4sin2 OwM)^' - (r|)™ - (r|), + «-'^(r|)^+" + "+'^(4)^-"] . (47b) 

Using the numbers from Eqs. ( [46a| through d) and Table 2, the IB contribution is apparently 
not negligible in either proton or neutron weak radii, and is indistinguishable from the 
strangeness radius contribution. The fact that both of the leading terms in Eq. (47a) are 
suppressed, i.e. that the weak neutral charge radius of the proton is naturally small, makes 
the forward angle proton asymmetry an attractive place to look for strangeness content. 
However, the relatively large size of the isospin breaking correction in this model, (^^-d- (^■f'^y+^ 
is more than 2% of the electromagnetic charge radius) shows that isospin corrections may 
prove to be a nontrivial effect which should be carefully taken into account. The shift in our 
calculated value of {r%)^^ due to mechanical IB is approximately +0.5% for both proton 
and neutron. For the proton, this shift is in fact +13% of the (very small) uncorrected 
Born approximation value .5(1 — 4 sin^ O^)- Even if our model result is viewed merely as an 
indication of the uncertainty introduced by isospin breaking, an extraction of {r'^)s by this 
means will be difficult at or below a level of ~ 2% of R^. This is equally relevant to the 
case of elastic '^He{e, e') experiments, where the isoscalar radius is measured. Here again, 
the isoscalar IB radius contributions add constructively (as would any intrinsic nucleon 
strangeness radii, of course). Note that the IB corrections due to the nuclear structure of 
'^He, i.e. without the intrinsic nucleon IB breaking, were recently calculated in Ref. |r^| and 
were found to be comparatively small. 



B. Electromagnetic Corrections 

As the first step in this analysis we take the lowest-lying baryon resonances with all the 
quantum numbers equal to those of the nucleon, but with "wrong" isospin as the leading 



14 



source of isospin admixtures. A number of such resonances can be found in e.g. Table 1.1 of 
Ref. all at the TV = 2 level in the harmonic oscillator model. Since the Coulomb potential 
does not carry orbital or spin angular momentum, each of these is separately conserved. That 
immediately eliminates all L = 2 states from this consideration. We are left with only one 
multiplet of L = resonances, the A(1550MeV), with the required properties^. 

The mass of A(1550) is anomalously low, i.e., at about the same mass as the M = 1 
negative-parity resonances. Irrespective of its present poorly confirmed experimental status, 
this is the only J\f = 2, SU(6) quark model state with the necessary quantum numbers; 
hence we will use it here. Its isospin-symmetric wave function is 



^10 : |A(1550)) 



^20 



1 



73 V TT 

V3 \~ 



exp 



2i?2 



2(p ■ A) exp 



2i?2 



{uud + udu + duu) 

v3 



{udd + dud + ddu) 



(48a) 
(48b) 
(48c) 
(48d) 
(48e) 



The next admixed terms occur at the M = A level, which ought to be suppressed due 
to the twice-as-large energy denominator. After some straightforward algebra we find the 
admixture coefficients from Eq. (|37D are 



-a 



3RAEV37T 



4.1 X 10" 



1 



--e„ = -2.0 X 10" 



(49a) 
(49b) 



where AE = En — -E'a(i55o) = —612 MeV is the experimental mass difference. 

At this stage we must say a few words about the choice of free parameters in this 
model. The constituent quark mass is fixed by fitting the nucleon magnetic moments, and 
is then entirely consistent with the mass of the nucleon. The oscillator frequency can be 
fixed in several ways, the two best known ones being: (a) using the observed spectrum of the 
nucleon resonances |T^, or (b) using the observed nucleon EM radius 0J^. The former leads 
to reasonable values of u only in the full-blown version of the model with strong hyperfine 
interaction and it underestimates the charge radius of the nucleon. The latter underestimates 
the energy gap between the ground and the excited states. While this discrepancy may be 
a serious cause of concern in attempts to make this model as realistic as possible, it is only 
mildly worrying in our case. Here, we are interested in a first estimate of isospin-breaking 
effects on electroweak form factors. We have therefore taken the former approach, and set 



In fact, this resonance does not appear in the latest Particle Data Group tables [16|, as it has 



been observed only by a single group |2C] 
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uo ^ 300MeV, which leads to (r^) = = (0.62 fm)^, an underestimate of the experiment by 
about 30%. This does correctly yield the observed A (1550) — mass difference, however. 
As one can see from Eqs. ( [49a| -b), an attempt to fix this up would somewhat decrease the 
admixture coefficients. In this sense, we are taking a conservative point of view, and our 
results might be viewed as a rough upper limit of the EM-induced isospin mixing effects. 

The appropriate EM transition matrix elements are all isoscalar, i.e., they are identical 
for the charged and neutral members of the multiplet 

(A|p(q)| N) = (A(1550)|f]e,exp(zq.ri)| N) 

(50a) 

\ UIILUJ I 

(A T |-^S^GM(q)| N]) = (A(1550) T | E -^^?e.exp(zq • r,)| N ]) 

exp ( —^1 (50b) 





72m V 3 V mu j \ Qmuj 
1 



(A T \jt'^'\ci)\ iV T) = (A(1550) T lE^^-rr exp(zq-r,)| N T) 

i=i ^ 

exp ( ) . (50c) 




\6muJ J 

None of these matrix elements survive the taking of the threshold limit. Hence, these 
corrections do not renormalize the static nucleon electroweak couplings. The only exception 
is the neutron anomalous magnetic moment which is renormalized by a small neutron EM 
mass shift (see below). The corrections to the effective radii are shown in Table III. The 
main conclusion of this section is that all of the explicit electromagnetic corrections due to 
admixtures with the A(1550) are extremely small, due to the smallness of a, coupled with 
the relatively small nuclear state overlaps, and relatively large energy denominator. 

In order to assess the reliability of our model, and that of our other assumptions, we 
apply them to a well-understood case of isospin breaking in nucleons: most notably the 
nucleon mass difference AM^ = Mp — Mn — —1.3 MeV. We calculate this mass difference 
with all free parameters in the model already fixed, compare it with experiment and then 
discuss other possible approaches (schemes) to fixing the free parameters. We find 

AMn = m^-md + 6Ec (51) 

where SEc = Eq — E^ = f (^rj^iT^j) is the difference of the Coulomb energies of the proton 
and the neutron. 

Using our Gaussian spatial wave functions, with R = 0.62 fm, the Coulomb energy 6Ec = 
■f^\l^ — 0.6 MeV. This reproduces the correct (negative) sign and order of magnitude (a 
few MeV) of the mass difference, but quantitatively our computed value of -3 MeV is more 
than two times too large. Our assumptions are also well known to lead to other unrealistic 
features of this model, such as Gaussian form factors. One can in principle increase the 



small Coulomb energy by turning on the strong hyperfine interaction between quarks, |14 
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but that would take us far outside of the intended scope of this article. We could also just 
reverse the procedure and fix the constituent quark mass difference from this calculation 
and the observed nucleon mass difference, but that would open new questions with regard 
to the applicability of those quark masses in the light meson sector. As far as the IB effects 
in the nucleon are concerned, that procedure would decrease their size as compared with our 
"direct" procedure. In that sense, as discussed earlier, we are again taking a conservative 
point of view. The above discussion ought to be an indication of the expected range of 
validity of this model: we may expect to have calculated the correct sign and order of 
magnitude of the effects of interest. We leave the investigation of next-order effects such as 
strong hyperfine interactions, and more realistic nucleon models as a task for the future. 



C. Discussion 



We would like to make a few comments with regard to the place of this investigation 
in the big picture of all "radiative corrections" that have to be evaluated. Our corrections 
form only one separately gauge-invariant class of purely EM radiative corrections. We 
have completely ignored EM quark vertex corrections and the associated self-energies and 
Bremsstrahlung effects. These form another separately gauge-invariant class of corrections 
that must be dealt with. If one assumes elementary i.e. Dirac constituent quarks, then one 
sub-part of such corrections, the anomalous magnetic moment of the quarks, are finite and 
gauge-invariant by themselves and are given by Schwinger's QED formula. This correction is 
manifestly of the same order of magnitude as the above-found IB nucleon-structure-induced 
effects. It is not clear to us whether they, or indeed any of the additional classes of graphs we 
have neglected, ought to be treated as intrinsic IB effects, or as pure radiative corrections. 
Furthermore, even if we were to include them into the former, the calculation would not be 
as straightforward as in free space, due to the strong interaction ("binding") effects. Similar 
considerations hold for the EM corrections to the axial current vertex. 

On the formal side, the present model is subject to the criticism that chiral symmetry 
is absent from it. If we assume that confinement and chiral symmetry breaking are two 
essentially independent phenomena, in accord with lattice QCD, then the situation can be 
remedied. The chiral symmetry and its spontaneous breakdown, which endows constituent 
quarks with their mass, can be installed into our model following e.g. Ref. Q], or Ref. pT] . 



These are, of course, just two of many chirally symmetric models for the constituent quark 
sector. Such chiral models allow a consistent evaluation of new kinds of diagrams such as 
the EM correction to the constituent quark axial current vertex shown in Fig. |^. Diagrams 
of this sort create a constituent quark isoscalar axial current coupling constant g^Qy which 
leads to an isoscalar nucleon axial current coupling constant gj^ ^ according to 

S S 3 S S 

^ = (iV t I^S^IAT t) = ^(A^ t T) = ^ . (52) 

2 2 2 2 

Manifestly, the result will depend on the model of the constituent quarks that one adopts. 
Should the vanishing of this class of IB corrections persist even at that next level of approx- 
imation, then we would be allowed the straightforward interpretation of a nonzero isoscalar 
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nucleoli axial coupling result, e.g. that of Ref. Q, in terms of strange quark contributions. 

V. CONCLUSIONS 

In summary, we have evaluated the IB admixture effects on the nucleon electroweak 
form factors in the nonrelativistic constituent quark model. We used the harmonic oscillator 
Hamiltonian as the confining interaction. The results were classified according to a model- 
independent formalism developed in Sect. II. 

If we are ultimately interested in either extracting strange nucleon form factors from 
electroweak data, or performing precision moderate energy Standard Model tests, we must 
first eliminate corrections such as the isospin breaking effects discussed here. These are 
certainly more "conventional" than the strangeness content, although also interesting and 
potentially important in their own right. A direct experimental measure of them, via e.g. 
Eq. (|T^, would be clearly valuable as well. Strange quark form factor estimates range con- 
siderably IQ, but typically yield from a small fraction of a percent up to > 10% corrections to 
electroweak structure. A priori, IB effects also modify electroweak structure of the nucleon, 
but only at a level of 0{a). The issue we have addressed here is effectively one of finding 
the coefficient in front of a in these corrections. If one based one's judgement solely on the 
observed nucleon mass splitting, then one would naively expect the correction to be very 
small, i.e. 0{AMn/Mn) ~ 0.1%, and not 1%. 

Our results show that the u-d mass difference can yield noticeably larger modifications 
to electroweak form factors than naively expected from the nucleon mass splitting. Never- 
theless, most of the corrections, with one possible exception, appear to be safely below the 
levels of what might be considered "interesting" strangeness content. Indeed, our results 
agree with the original expectation of 0{a) corrections, as well as with our assumption that 
the EM-mechanical interference (cross-) terms may be neglected. In the isoscalar charge 
radius, we find corrections of several percent. The main source of this IB is simply the 
shift of the nucleon's center of mass. It is in this regard that the present model is very well 
suited to the task at hand, in contrast to e.g. the bag models. We have, however, chosen 
the model parameters to maximize these IB effects, so within the subclass of IB corrections 
considered here, our numbers should represent an upper limit. E.g., choosing Arriq to re- 
produce the observed Mp — Mn splitting in our model would halve the predicted mechanical 
isospin corrections to the charge radius. Since IB breaking can, in principle, interfere with 
the extraction of a relatively modest strangeness-carrying nucleon electric form factor from 
lepton-nucleon scattering experiments, our results indicate the need for careful checking of 
these results in other nucleon models. 



^ There are no odd parity admixture corrections to leading order in the experiment since that 
is a neutrino scattering experiment. 
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FIGURES 



FIG. 1. Feynman diagrams describing the EM isospin-breaking contribution to the polar- and 

axial- vector currents. The solid line denotes the quarks, wavy solid line is the photon and wiggly 
solid line is the neutral intermediate vector boson Z^. The shaded blob together with the three 
solid lines and one double solid line leading to it denotes the nucleon wave function. 

FIG. 2. An example of an 0{a) EM correction to the quark EM and axial current vertices 
of the type that may renormalize the nucleon magnetic moment and the axial coupling constant. 
These diagrams are not evaluated in this paper. 
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TABLES 



TABLE L The shift of nucleon static properties due only to mechanical IB corrections in our 
model. Here M„ = 2mu + ma; = 2mii + and x is defined by x = sin^ 9y/- Numerical values 



used are given before Eq. 19 in the text. 
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TABLE II. Theoretical nucleon radii (r^) with mechanical IB corrections, as well as the relative 
change. Notation is defined in the text, after Eqs. (44a - d), except x = siv?Ow- The isospin 
symmetric nucleon radius in each case is (except (r|)™ = 0) where we use = .62 fm. 
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TABLE III. The shift of nucleon radii 5(r^) due to purely electromagnetic IB corrections, as 



calculated in section 4B. Here £p = 4.1 x 10 e„ = —1/2 Sp. 
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